Defining a robust biological prior from Pathway Analysis to drive Network Inference

Inferring genetic networks from gene expression data is one of the most challenging work in the post-genomic era, partly due to the vast space of possible networks and the relatively small amount of data available. In this field, Gaussian Graphical Model (GGM) provides a convenient framework for the discovery of biological networks. In this paper, we propose an original approach for inferring gene regulation networks using a robust biological prior on their structure in order to limit the set of candidate networks. Pathways, that represent biological knowledge on the regulatory networks, will be used as an informative prior knowledge to drive Network Inference. This approach is based on the selection of a relevant set of genes, called the "molecular signature", associated with a condition of interest (for instance, the genes involved in disease development). In this context, differential expression analysis is a well established strategy. However outcome signatures are often not consistent and show little overlap between studies. Thus, we will dedicate the first part of our work to the improvement of the standard process of biomarker identification to guarantee the robustness and reproducibility of the molecular signature. Our approach enables to compare the networks inferred between two conditions of interest (for instance case and control networks) and help along the biological interpretation of results. Thus it allows to identify differential regulations that occur in these conditions. We illustrate the proposed approach by applying our method to a study of breast cancer's response to treatment

A meta-analysis of randomized double-blind clinical trials in CMT1A to assess the change from baseline in CMTNS and ONLS scales after one year of treatment

International audienceCMT1A is the most common inherited peripheral neuropathy. There is currently no approved treatment. We performed a meta-analysis including four randomized, double-blind, Placebo-controlled clinical trials to assess the disease progression after one year under Placebo, Ascorbic Acid (AA) or PXT3003, a combination of three repurposed drugs. We observed a weak deterioration in patients under Placebo, well below the reported natural disease progression. Patients treated with AA were stable after one year but not significantly different from Placebo. Patients undergoing PXT3003 treatment showed an improvement in CMTNS and ONLS, statistically significant versus Placebo and potentially precursory of a meaningful change in the disease course

Unraveling the relationships between alpha- and beta-adrenergic modulation and the risk of heart failure

Background: The effects of α and ß adrenergic receptor modulation on the risk of developing heart failure (HF) remains uncertain due to a lack of randomized controlled trials. This study aimed to estimate the effects of α and ß adrenergic receptors modulation on the risk of HF and to provide proof of principle for genetic target validation studies in HF. Methods: Genetic variants within the cis regions encoding the adrenergic receptors α1A, α2B, ß1, and ß2 associated with blood pressure in a 757,601-participant genome-wide association study (GWAS) were selected as instruments to perform a drug target Mendelian randomization study. Effects of these variants on HF risk were derived from the HERMES GWAS (542,362 controls; 40,805 HF cases). Results: Lower α1A or ß1 activity was associated with reduced HF risk: odds ratio (OR) 0.83 (95% CI 0.74–0.93, P = 0.001) and 0.95 (95% CI 0.93–0.97, P = 8 × 10−6). Conversely, lower α2B activity was associated with increased HF risk: OR 1.09 (95% CI 1.05–1.12, P = 3 × 10−7). No evidence of an effect of lower ß2 activity on HF risk was found: OR 0.99 (95% CI 0.92–1.07, P = 0.95). Complementary analyses showed that these effects were consistent with those on left ventricular dimensions and acted independently of any potential effect on coronary artery disease. Conclusions: This study provides genetic evidence that α1A or ß1 receptor inhibition will likely decrease HF risk, while lower α2B activity may increase this risk. Genetic variant analysis can assist with drug development for HF prevention

A Reliable Method for the Selection of Exploitable Melanoma Archival Paraffin Embedded Tissues for Transcript Biomarker Profiling

The source tissue for biomarkers mRNA expression profiling of tumors has traditionally been fresh-frozen tissue. The adaptation of formalin-fixed, paraffin-embedded (FFPE) tissues for routine mRNA profiling would however be invaluable in view of their abundance and the clinical information related to them. However, their use in the clinic remains a challenge due to the poor quality of RNA extracted from such tissues. Here, we developed a method for the selection of melanoma archival paraffin-embedded tissues that can be reliably used for transcript biomarker profiling. For that, we used qRT-PCR to conduct a comparative study in matched pairs of frozen and FFPE melanoma tissues of the expression of 25 genes involved in angiogenesis/tumor invasion and 15 housekeeping genes. A classification method was developed that can select the samples with a good frozen/FFPE correlation and identify those that should be discarded on the basis of paraffin data for four reference genes only. We propose therefore a simple and inexpensive assay which improves reliability of mRNA profiling in FFPE samples by allowing the identification and analysis of “good” samples only. This assay which can be extended to other genes would however need validation at the clinical level and on independent tumor series

Clinically Relevant Characterization of Lung Adenocarcinoma Subtypes Based on Cellular Pathways: An International Validation Study

Lung adenocarcinoma (AD) represents a predominant type of lung cancer demonstrating significant morphologic and molecular heterogeneity. We sought to understand this heterogeneity by utilizing gene expression analyses of 432 AD samples and examining associations between 27 known cancer-related pathways and the AD subtype, clinical characteristics and patient survival. Unsupervised clustering of AD and gene expression enrichment analysis reveals that cell proliferation is the most important pathway separating tumors into subgroups. Further, AD with increased cell proliferation demonstrate significantly poorer outcome and an increased solid AD subtype component. Additionally, we find that tumors with any solid component have decreased survival as compared to tumors without a solid component. These results lead to the potential to use a relatively simple pathological examination of a tumor in order to determine its aggressiveness and the patient's prognosis. Additional results suggest the ability to use a similar approach to determine a patient's sensitivity to targeted treatment. We then demonstrated the consistency of these findings using two independent AD cohorts from Asia (N = 87) and Europe (N = 89) using the identical analytic procedures

Should We Abandon the t-Test in the Analysis of Gene Expression Microarray Data: A Comparison of Variance Modeling Strategies

High-throughput post-genomic studies are now routinely and promisingly investigated in biological and biomedical research. The main statistical approach to select genes differentially expressed between two groups is to apply a t-test, which is subject of criticism in the literature. Numerous alternatives have been developed based on different and innovative variance modeling strategies. However, a critical issue is that selecting a different test usually leads to a different gene list. In this context and given the current tendency to apply the t-test, identifying the most efficient approach in practice remains crucial. To provide elements to answer, we conduct a comparison of eight tests representative of variance modeling strategies in gene expression data: Welch's t-test, ANOVA [1], Wilcoxon's test, SAM [2], RVM [3], limma [4], VarMixt [5] and SMVar [6]. Our comparison process relies on four steps (gene list analysis, simulations, spike-in data and re-sampling) to formulate comprehensive and robust conclusions about test performance, in terms of statistical power, false-positive rate, execution time and ease of use. Our results raise concerns about the ability of some methods to control the expected number of false positives at a desirable level. Besides, two tests (limma and VarMixt) show significant improvement compared to the t-test, in particular to deal with small sample sizes. In addition limma presents several practical advantages, so we advocate its application to analyze gene expression data

Accounting for Population Stratification in Practice: A Comparison of the Main Strategies Dedicated to Genome-Wide Association Studies

Genome-Wide Association Studies are powerful tools to detect genetic variants associated with diseases. Their results have, however, been questioned, in part because of the bias induced by population stratification. This is a consequence of systematic differences in allele frequencies due to the difference in sample ancestries that can lead to both false positive or false negative findings. Many strategies are available to account for stratification but their performances differ, for instance according to the type of population structure, the disease susceptibility locus minor allele frequency, the degree of sampling imbalanced, or the sample size. We focus on the type of population structure and propose a comparison of the most commonly used methods to deal with stratification that are the Genomic Control, Principal Component based methods such as implemented in Eigenstrat, adjusted Regressions and Meta-Analyses strategies. Our assessment of the methods is based on a large simulation study, involving several scenarios corresponding to many types of population structures. We focused on both false positive rate and power to determine which methods perform the best. Our analysis showed that if there is no population structure, none of the tests led to a bias nor decreased the power except for the Meta-Analyses. When the population is stratified, adjusted Logistic Regressions and Eigenstrat are the best solutions to account for stratification even though only the Logistic Regressions are able to constantly maintain correct false positive rates. This study provides more details about these methods. Their advantages and limitations in different stratification scenarios are highlighted in order to propose practical guidelines to account for population stratification in Genome-Wide Association Studies